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Abstract: We present hybrid block of eight integrators which are of uniform order nine through interpolation 
and collocation procedures. The block method is implemented on a stiff, physical and life problems in which 
the results obtained compares favourably with the existing method. The properties of the hybrid block 
integrators is fully investigated and confirmed to be computationally stable with the numerical experiments 
tested. 


1.0 Introduction 

Among the most important mathematical tools used in producing models in the physical 
sciences, Biological sciences and Engineering are differential equations. But most of these 
differential equations do not posses closed form or finite solutions. Even if they posses 
closed form solutions we do not know the method of getting them. 

In many real-life situations, the differential equation that models the problem is too 
complicated to solve exactly. Hence there is need to develop an accurate algorithms for 
obtaining an equivalent approximating solution to the original problems. Most recent 
researchers has developed some block methods to cater for this class of problems 
y=f@y),  y(a-p (1) 
Among such researchers are [1], [2], [3], [4], [5], [6] and [7] to mention a few. 

In our method, a block of eight integrators were proposed at step length of four which are 
of uniform order nine and also all the Discrete schemes in our block came from a single 
continuous formula. 


2.0 | Development of the method 
Our objective is to derive a block of eight integrators of the form 


Y(Xn+v) =X0 Yn + h Bohr + EPT + [fusi EME + b2fn+2 + Bft + B3fn+3 


+ furit Ps] (2) 


3 3 5 7 
where Xn: Y 2 2,1,2,2, 2,357, 4. 


We proceed by seeking an approximate solution of the form 
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mt+t-1 


yx) = 5, wx! (3) 
j70 

y= jay xt = f(xy) (4) 
j=l 


where m and t are the number of collocation and interpolation points used in the method. 
Specifically for this method m = 9, t = 1 and the degree of the polynomial is m + t — 1. 


Equation (3) is interpolated at x = x, and (4) is collocated at x = x,,,, v = 


0,2 , 1, =, 2, 5 3, 3 and 4 which leads to the following non linear system of equations of 


the form 
mt+t-1 
P(x) = 2. oc; xj) 
j-0 
mt+t-1 mt+t-1 
PO)- $ Foa =h DB Of = fey) (5) 
j=1 j=0 


When using Maple 17 (Mathematical software) to determine the unknown parameters c; 


and f; in (5), we obtain the following. 
Xo= Yn 
T8027 p 71277 2732 6989 1369 161 
75601. + 2268012 1080}? 5400h4 3240h” 1890h6 
aD 4 1? 
7560h7 8505A? 
"T 4096(2520) , 8192(8658) , , 2048(1745) , 8192(1316) . , 4096(575) |, 
a 405405h 12162157? 57915h3 289575h4 173745h5 
16384(73) , 4096(35) , ^ 65536 


^. 405405h$ 405405h6  3648645h8 





9 
































7560 28494 , 42783 33713 . 7605 2(989) 69 
B= \>aeq) i ea Gea l6 + 7 -— 18 
180h 270h? 360h3 450h4 270h° 315h6 90h6 
di r| 
405h8 
pe Pes jp. 4(10338) 3 _ 1(16867) , 4(1425) . | 2(6805) , 16(463) , | 2(67) ja 
2 135A 405A? 135h3 675h* 405h5 945h® 135h® 
64 
-iugs"l 
1215/8 
2 Li ;,1(016137) , 1(13785) , | 124463) . 1(6115) , 1(867) , 
* 60h 180h2 120h3 300h* 180A5 1055 
1(65) , 
LL r| 
60h6 ^s 
N POST?) 2 _ 8(6606) ,  2(1659) , 8(1522) s , 4(789) js 32(58) 14 4(63) ja 
p 315A (945h2 45h3  225h^. 135h5 315h6 315h6 


128 i 
28359 





























A NEW NINETH ORDER HYBRID BLOCK INTEGRATOR 3/8 
_ e 2 , 1(11178)  _ 1(20033) , _ 1(18821) , 1(5017), | 2(761) , 
3 324h 486h2 648h? 810h4 486h5 567h® 
BC: J 
162h6 729h8 
jos 2(1080),, 4(4842),,  1(1257),, 4(1202), 2(655), 16(51),, 269) 1g 
Z 1155h 3465h2 165h3 825h* 495h5 1155/5 1155h6 
64 
-——Ó 
103958 
L[ 10215),  1(8541),, 1(3921),  1(15203), 1(4215) , | 1(671) j, 
i 14040h 14040h2 4680h? 23400n^ 14040A5 8190h® 
1(57) y 4 
= r| 
4680h® 5265h8 
(6) 
where | = (x — x4) 
y(x) 2«o Yn +h Bohr + Bsf 43 + Bifnti + Bsf 43 + Bofn+2 + Bsf 5 + Bsfn+3 
4 "+3 Zz ntz 2 END 
+ Bf ar Bafnsa| (7) 
Pe 


Equation (6) is substituted in equation (7) to obtain our continuous formula. Also 


3 3 5 


evaluating (7) at x 2x,,, V = n 1,5 id ^ 3 E and 4 to yield eight integrators to form 


our hybrid block methods as follows: 
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00010 0 0 ol»2| 10 00 0 0 0 0 1l| Ya 
000010 0 0|X,s| |0 0 0 0 0 O 0 1[||X,.3 
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1925 89600 22400 ~ 179200 1600 17920 
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+n) 18243225 — 14175 42525 ~ 4725 2025 368550 
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729729 290304 217728 ~ 193536 36288 2612736 
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(8) 


3.0 Basic Analysis of the proposed hybrid block method 


3.1 Order and Error constants ofthe newly hybrid block integrators according to [8], 
given a linear multi-step method of the form 


k k 
>, yma = A9 By GO; (9) 
j-0 j-0 
We associate the linear difference operator 
k 
£[y GO] = > [sy yGe + th) — hgy (x + jh] (10) 
j=0 


where y(x) is an arbitrary function, continuously differentiable on [a, b]. the test function 
Expanding y(x + ih) and y' (x + jh) as Taylor series about x and collecting like terms in (10) 
gives 

L[y(x), h] = coy (x) + cihy'® += + co hy (x) (11) 
Definition 3.1.1 
A linear multistep method of the form 


mt+t-1 mt+t-1 
2, Gy =h > By Casi 
j=0 j-0 
is said to be of order p, if from (11), we have cp = c4 = c; = + Cp = 0 but C,,, # 0 and 


Cp+1 is called error constant. [8] 


Definition 3.1.2 
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A LMM of (9) is said to satisfy the root conditions if all of the roots of the first 
characteristics polynomial have modulus less than or equal to unity and those of modulus 
unity are simple. The method (9) is said to be zero-stable if it satisfies the root condition. [8] 


Table 1 Order and error constants of integrators (8) 





























Schemes Order Error Constants 
V3 9 1627529 
nat — ——— — 
93952409600 
ya 9 1027 
59719680 
Ypa? 9 6439 
> PRESA NER 
367001600 
Yna2 9 2257 
130636800 
ys 9 187975 
nt, "TI WEIN X M 
10701766656 
Yn«3 9 7 
409600 
y 7 9 140581 
ng ee 
7644119040 
Vn+4 9 23 
2041200 

















4.0 Consistency and Stability of the block schemes 
The hybrid block integrators (8) is said to be zero stable, if the roots Z, ,S = 1,2,.. k of 
the first characteristic polynomial p(z) = det A(49 — AQ) satisfies |z,| € 1 has 


multiplicity not exceeding the order of the differential equation. 


10000000 0000000 1 
0 1 0 0 00 9 9 0000000 1 
00100000 0000000 1 

s od OL 50:20-22:09: 29559: 9 |, 97 30730" 20: 0: 50: NII. 

"2-Mlooo0o010200 00000020210| 9 
00000100 0000000 1 
00000010 00000001 
00000001 00000001 

AA-1) =0 


A4 = 0,45 = 0,45 = 0,44 = 4, Às = 0, Ag = 0,4; = 0, Ag = 1 
From definition 3.1.2, the newly hybrid block method (8) is zero stable and also 
consistent since the order of all the integrators in (8) are 9 > 1 
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5.0 | Numerical experiments 
The hybrid block method are demonstrated on stiff, life and physical problems 
Example 5.1 

Qe ey 

y(0) =0,h=0.1 

Analytic solution: iua 

e?*41 

Example 5.2 (SIR Model) 

The SIR model is an epidemiological model that computes the theoretical number of 
people infected with a contagious illness in a closed population over time. The name of 
this class of models derives from the fact that they involve coupled equations relating the 
number of susceptible people S(t), number of people infected /(t), and the number of 
people who have recovered R(t),. This is a good and simple model for many infectious 
diseases including measles, mumps and rubella (Rosser 1967). It is given by the following 
three coupled equations. 


dS 
pr ee 


or Tanhx 





dt 
dl 
di^ ul — yl + BIS (12) 
dM R t yl 
where yg, y and f are positive parameters. Define y to be, 

y=S+I+R (13) 
When solutions in (12) are substituted in (13) we have 
y = u(1 — y)t then 

y'=u-y) y@=p (14) 
Taking uw = 0.5 and attaching an initial condition y(0) = 0.5 (for a particular closed 
population), we obtain 

y —-05(1—-y) yO)=0.5 (15) 


Exact solution: y(t) = 1 — 0.5e 9t 


Table 2 Performance of Hybrid block method (8) on problem 5.1 






































x Exact solution Computed solution Error of [9] New Error 
0.1 0.099667994624956 0.099667994698116 1.2312 E(-09) 7.31 E(-11) 
0.2 0.197375320224904 0.197375320296238 1.2343 E(-09) 7.1334 E(-11) 
0.3 0.291312612451591 0.291312612518669 1.2814 E(-09) 6.7078 E(-11) 
0.4 0.379948962255225 0.379948962292507 2.2283336 E(-09) 3.7282 E (C11) 
0.5 0.462117157260010 0.462117157294015 3.3949 E(-09) 3.4005 E (-11) 
0.6 0.537049566998035 0.537049567020542 3.5428 E(-09) 2.2507 E (-11) 
0.7 0.604367777117164 0.604367777137633 3.8244 E(-09) 2.0469 E (-11) 
0.8 0.664036770267849 0.664036770297224 3.8293 E(-09) 2.9375 E (-11) 
0.9 0.716297870199024 0.716297870224885 1.3275 E(-09) 2.5861 E (-11) 
1.0 0.761594155955765 0.761594155986542 1.4427 E(-09) 3.0777 E (-11) 
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Figure 1: Error graph of solution of example 5.1 


—$-— Error of [9] 


== New Error 


Table 3 Performance of Hybrid block method (8) on problem 5.2 
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t | Exact solution Computed solution Error of [3] New Err 

0.1 | 0.524385287749643 | 0.524385287749643 | 5.57443 E(-12) | 0.00 

0.2 | 0.547581290982020 | 0.547581290982018 | 3.946177 E(-12) | 2.00 E(-15) 
0.3 | 0.569646011787471 | 0.569646011787471 | 8.183232 E(-11) | 0.00 

0.4 | 0.590634623461009 | 0.590634623461001 | 3.436118 E(-11) | 8.00 E (-15) 
0.5 | 0.610599608464298 | 0.610599608464286 | 1.92974 E(-10) 1.20 E (-14) 
0.6 | 0.629590889659141 | 0.629590889659125 | 1.87904 E(-10) 1.60 E (-14) 
0.7 | 0.647655955140644 | 0.647655955140626 | 1.776835 E(-10) | 1.80 E (-14) 
0.8 | 0.664839976982180 | 0.664839976982157 | 1.724676 E(-10) | 2.30 E (-14) 
0.9 | 0.681185924189114 | 0.681185924189090 | 1.847545 E(- 2.40 E (-14) 

10) 
1.0 | 0.696734670143684 | 0.696734670143655 | 3.00577 E(-10) 2.90 E (-14) 














—$-— Error of [3] 


== New Error 


Figure 2: Error graph of solution of example 5.2 
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Discussion of Results 
We observed that the newly block integrator (8) performed excellently well with the 
two problems tested with the schemes. This shows that our method is good and can be 
used to solve accurately any model of the form y' = f(x,y), y(a) = p. (see tables 2,3 
and figures 1,2) 


Conclusion 
We want to conclude that the newly block integrator (8) is of uniform order 9, zero 
stable, consistent and self starting. Also result obtained from the method converges 
more than the existing method. (see figures 1 and 2) 
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